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Abstract 

This  paper  explores  dwell  time  constraints  on  switched  systems  with  mul¬ 
tiple,  possibly  disparate  invariant  limit  sets.  We  show  that,  under  suitable 
conditions,  trajectories  globally  converge  to  a  superset  of  the  limit  sets  and 
then  remain  in  a  second,  larger  superset.  We  show  the  effectiveness  of  the 
dwell-time  conditions  by  using  examples  of  switching  limit  cycles  commonly 
found  in  robotic  locomotion  and  flapping  flight. 
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1.  Introduction 

Bifurcations  have  been  of  interest  to  dynamical  systems  theory  for  decades. 
However,  most  control  strategies  view  such  behavior  as  damaging  and  try  to 
mitigate  it  [1].  Relatively  less  work  actively  inserts  bifurcations  as  part  of  a 
control  strategy.  One  example  is  using  a  classic  Hopf  bifurcation  for  mode¬ 
switching  between  flapping  and  gliding  flight  in  micro-aerial  vehicles  [2], 
The  authors  consider  a  supercritical  Andronov-Hopf  bifurcation  model  of 
x  =  («;  v): 


/—A/p2  (u2  +  v2  —  p2a )  u  —  oj(t)v\ 
y  (jj{t)u  —  A/p2  (u2  +  v2  —  p2a)  v  J 
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Figure  1:  Schematic  of  mode  switching  with  non-equilibrium  limit  sets. 


with  cr  =  1.  For  a  positive  rate  of  convergence  A  >  0,  it  can  be  easily  shown 
that  any  initial  trajectory  («;  v)  /  0  exponentially  converges  to  a  circle 
of  the  radius  p  rotating  at  the  time- varying  frequency  uj(t)  with  bounded 
uj{t).  If  a  <  0,  bifurcation  occurs  and  the  system  globally  converges  to  the 
origin,  which  is  useful  for  fast  inhibition  of  oscillation.  Fast  inhibition  and 
synchronization  of  oscillators  are  key  properties  for  many  neurobiologically- 
inspired  control  schemes. 

Another  possible  application  is  walking  robots.  Fig.  1  shows  a  hypotheti¬ 
cal  switching  pattern  for  a  walking  robot  application  utilizing  central  pattern 
generation.  A  guidance/navigation  engineer  may  design  limit  cycle  subsys¬ 
tems  for  walking  and  jumping  modes  (shown  as  a  Hopf  oscillator  and  a  Van 
der  Pol  oscillator),  while  utilizing  steady-state  control  strategies  for  static 
balancing  or  tasks  requiring  fine  motor  control. 

Mode-switching  also  implicates  a  large  body  of  literature  on  switched 
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systems  [3].  Most  work  on  stability  of  switched  systems  assumes  that  all 
subsystems  have  a  common  equilibrium  point.  [4,  5,  6]  consider  weak  Lya¬ 
punov  functions  in  the  style  of  LaSalle  for  a  common  equilibrium.  [7]  consid¬ 
ers  equilibrium  location  changes,  but  holds  the  vector  field  constant.  They 
connect  the  result  to  averaging  theory.  [8]  considers  practical  stability  of 
affine  systems  with  multiple  distinct  equilibria.  Alpcan  and  Ba§ar  investi¬ 
gated  dwell  time  criteria  for  nonlinear  globally  exponentially  stable  subsys¬ 
tems  which  could  have  differing  equilibria  [9].  Such  systems  have  no  single 
globally  attractive  equilibrium  point.  The  authors  of  [9]  reported  an  explicit 
construction  of  the  dwell  time  and  a  conservative  invariant  set.  This  paper 
is  inspired  by  that  work  and  is  a  generalization  of  it.  We  generalize  their 
result  to  switched  systems  where  each  subsystem  may  have  multiple  invari¬ 
ant  sets.  We  pursue  a  similar  dwell  time  strategy  in  order  to  provide  spatial 
bounds  for  the  switched  system.  The  resulting  construction  is  slightly  more 
complicated,  as  we  consider  V  in  order  to  isolate  the  invariant  sets  rather 
than  using  the  Lyapunov  function  alone. 

Systems  with  bifurcation  often  contain  multiple  cu-limit  sets  which  cannot 
be  globally  exponentially  stable.  Instead,  results  such  as  LaSalle’s  invariant 
set  theorem  [10]  allow  us  to  analyze  asymptotic  stability  of  this  larger  class 
of  systems.  LaSalle’s  theorem  and  much  of  the  switched  systems  literature 
are  both  Lyapunov-based,  and  we  will  make  use  of  Lyapunov  functions  to 
define  all  the  relevant  sets.  The  benefit  to  relying  on  Lyapunov  functions  is 
that  this  requires  no  special  structure  on  the  subsystems’  entire  vector  fields. 
The  tradeoff  is  that  we  fail  to  exploit  any  special  structure  the  subsystems 
may  possess  and  the  result  relies  on  being  able  to  find  suitable  Lyapunov 
functions. 

Section  2  provides  background  assumptions  and  definitions.  Section  3 
begins  by  reconsidering  existing  results.  Sections  3.3  through  3.5  present  two 
methods  to  accomplish  the  goal.  Choice  of  a  particular  method  will  depend 
on  specific  situations  and  design  constraints.  Section  4  shows  a  numerical 
example,  and  Section  5  provides  closing  comments. 

2.  Preliminaries  and  Definitions 

Consider  a  set  of  continuous-time  dynamical  systems  defined  by 

x  =  fp  (x) ,  (2) 
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where  x  G  Mn  and  p  G  V,  with  some  index  set  V  =  {pi,P2,  •••,Pmax}-  A 
piecewise  constant  switching  signal  a  :  [0,  oo)  — >  V  specifies  the  active  sub¬ 
system  at  each  time.  Assume,  for  ease  of  analysis,  that  fp  are  each  continuous 
with  continuous  first  partials.  Together,  (2),  the  index  set,  and  the  switching 
signal  define  a  switched  system. 

Only  some  systems  admit  stability  results  for  arbitrary  switching  signals, 
so  we  will  consider  a  constraint  on  how  quickly  the  switching  signal  can  make 
consecutive  switches. 

Definition  Consider  a  switched  system  with  switching  times  It 

is  said  to  have  dwell  time  r  if  U+i  —  U  >  r  Vi  G  N. 

Next,  we  review  and  introduce  some  important  subsets  of  Rn.  We  have 
not  yet  provided  strict  assumptions  on  Lyapunov- like  functions.  At  present, 
it  is  enough  to  assume  that  each  subsystem  has  a  (possibly  different)  C 1 
Lyapunov-like  function,  which  is  bounded  above  and  below  on  every  bounded 
subset  of  Mn.  Furthermore,  assume  that  each  is  radially  unbounded  (V^(x)  — >• 
oo  as  ||x||  — y  oo).  This  ensures  that  every  sublevel  set  describes  a  compact 
region.  We  assume  for  the  remainder  of  this  paper  that  the  minimum  value 
of  each  Vp  is  zero.  Define 


£p  =  {xGMn|Dp(x)  =  0}  (3) 

as  the  set  which  attains  the  minimum  value  of  Vp.  Let  k  be  a  positive  constant 
and  define 


Np(k)  =  {x  G  Mn|Dp(x)  <  «},  (4) 

a  closed  ^-neighborhood  of  Qp.  For  the  purposes  of  Theorem  1,  Afp  (k)  is 
connected,  but  it  is  not  necessarily  connected  in  the  remainder  of  the  paper 
(See  Figure  4).  Let  the  union  over  V  be 

A /"(«)  =  (JA/^k).  (5) 

pev 

Additionally,  we  define  a  superset,  A4(k),  in  a  series  of  steps  with 


ap  (k)  =  max  Vp  (x) , 

xGA/^k) 

and 

Mp  (k)  =  {x6rTp(x)<  ap  (k)}  . 


(6) 

(7) 
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Figure  2:  Qualitative  example  of  how  A f  and  M  are  built  for  a  switched  system  consisting 
of  two  subsystems,  each  with  a  single  equilibrium,  but  at  different  locations. 

Finally,  we  create  a  closed  union  of  closed  sublevel  sets, 

M(k)=\Jmp.  (8) 

pev 

Notice  that  the  dependence  on  k  carries  through  once  we  use  it  in  J\f  (n).  For 
the  purposes  of  Theorem  1,  M.  is  a  connected  superset  of  A/".  Theorem  4 
will  introduce  a  different  notion  which  is  not  necessarily  connected.  Figure  2 
provides  a  one-dimensional  example  to  help  visualize  these  sets. 

3.  Stability  Results 

We  first  restructure  the  result  in  [9]  slightly  to  add  clarity  and  to  better 
facilitate  the  generalization  presented  in  this  paper. 

3.1.  Unique  Equilibrium  Case 

Theorem  1  ([9]).  Consider  a  family  of  systems  defined  by  (2),  each  with  a 
single,  globally  exponentially  stable  equilibrium,  denoted  x*.  Suppose  that  the 
exponential  decay  rate  of  each  Lyapunov  function,  as  described  in 

Vp(x)  <  -eVp(x),  Vx  G  IT,  Vp  G  V,  (9) 

is  at  least  e  >  0. 

Furthermore,  given  a  positive  constant  k,  define  the  sets  as  in  Section  2 
and  assume  / i(k )  G  (l,oo)  such  that 
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Then,  for  every  switching  signal  with  dwell  time 


log  h(k) 

T  >  - > 

£ 

(i)  There  exists  a  time  T  such  that  x(T“)  G  A C(t-)  (k),  and 

(ii)  For  any  time  t  such  that  x(t)  G  Afa{t)  («),  x(t)  G  Ad  (/c)  /or  all  t  >  t. 

Remark  1.  Not  all  choices  of  Vp  will  give  a  finite  value  of  /i(k).  Typical 
polynomial  constructions  for  Lyapunov  functions  must  have  the  same  poly¬ 
nomial  order.  Scaling  or  stretching  Lyapunov  functions  may  be  useful,  but 
there  will  be  implications  for  both  the  spatial  parameter  /i  and  the  temporal 
parameter  e.  We  will  will  revisit  the  idea  of  scaling  briefly  in  Sections  3. 3-3. 4. 

Proof.  We  will  only  provide  a  sketch  of  the  relevant  features.  The  proof 
proceeds  in  two  parts: 

(i)  Consider  a  finite  time  interval  [f0,  T]  with  corresponding  switching  times 
ti  <  t2  <  •  •  •  <  tn<r ,  where  na  is  the  number  of  switches  inside  the 
interval.  Between  switches,  aft)  is  constant.  If  the  trajectory  enters 
A/r(t)  (not  just  A/"),  the  result  is  trivial.  Otherwise,  the  behavior  of 
x(t))  between  switches  satisfies  (9).  Denote  the  limit  from  the 
right/left  as  superscript  +/— ,  respectively.  Then, 

<  e-£(‘1+1-‘l)U(t?|(x(t+)).  (12) 

At  switches, 

^(t+)(x(A+))  <M«)K(tn(x(*r))  (13) 

holds.  We  can  iterate  on  i  to  obtain 

W(T-)(X(T-)  <  e(OogM«)/r)-e)(T-to)yff(to)(x(to)).  (14) 

Importantly,  under  the  dwell  time  condition,  (11),  this  implies  that  by 
taking  T  suitably  large,  we  can  make  Va(r-)fx.{T~))  arbitrarily  small. 
Thus,  x(T_)  G  A fa{T~)  (K)-  This  proof  is  existential,  not  constructive. 
We  cannot  calculate  a  particular  time  T  for  any  particular  problem. 


(ii)  The  second  part  of  the  proof  shows  that  after  a  switch  at  time  tt,  the 
dwell  time  is  sufficiently  large  to  force  the  trajectory  back  into  J\fa(t)  («) 
before  a  subsequent  switch  at  time  U+ Furthermore,  the  trajectory 
cannot  escape  Ad  (k)  in  that  interval. 

□ 

Details  are  available  in  [9].  The  proof  presented  above  is  different  from  [9] 
in  one  important  way  -  we  specify  that  the  trajectory  enters  JVa(t)  (K)  rather 
than  In  fact,  it  is  an  error  to  do  the  latter.  The  trajectory  may 

pass  through  _A/"(k)  \  (K)  and  then  switch  after  it  has  emerged,  which 

may  cause  it  to  exit  Ad  (n).  There  is  nothing  special  about  \J\fa(t)  («) 

while  cr(t)  remains  constant.  It  may  be  aesthetically  displeasing  to  have  the 
entry  set  change  in  time  along  with  a (t),  but  we  must  do  this,  because  the 
time-independent  formulation  is  false.  We  will  demonstrate  with  an  example. 

Example  This  example  is  a  slight  variation  of  Example  2  from  [9] .  Choose 

xp  =  Axp  +  bp,  (15) 

but  with 


A 


-1  -10  ' 
10  -1 


(16) 


We  are  able  to  use  the  same  Lyapunov  functions,  El(x)  :=  x\  +  (.xq  —  l)2, 
V2 (x)  :=  (xi  +  l)2  +  .h],  and  Vs(x)  :=  {%i  ~  l)2  +  One  can  check  that 
r  can  be  the  same  as  in  [9],  but  that  is  not  important  here.  Consider  the 
trajectory  shown  in  Figure  3.  The  red  circles  show  A/"(k);  the  black  circles 
show  Ad  (k).  Subsystem  1  is  active  at  the  start,  and  the  trajectory  is  shown 
in  blue.  A/"i  (k)  is  the  red  circle  centered  around  [  0  1  ]T.  The  trajectory 

passes  through  A3  (/-c),  which  is  the  red  circle  centered  around  [  1  0  ]T. 
This  is  entering  A /"(ft),  but  not  Afa(t)  (K)-  Now,  notice  that  we  could  have 
selected  initial  conditions  to  make  the  entry  into  A3  (k)  occur  at  arbitrarily 
large  time.  This  allows  us  to  place  a  switch  anywhere  along  the  trajectory, 
regardless  of  r  (this  is  a  “free  switch”  that  we  will  see  again  in  Section  3.5). 
In  this  case,  we  switched  to  subsystem  2  at  the  location  where  the  trajectory 
changes  to  green.  It  exits  Ad  A)  shortly  thereafter. 
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System  Trajectory  x(t) 


Figure  3:  Example  trajectory  demonstrating  the  need  for  time-dependent  Afa(t)  («)■ 


3.2.  Problem  Statement:  Switching  Systems  having  Multiple  Invariant  Sets 

Alpcan  and  Ba§ar  considered  subsystems,  each  having  a  globally  expo¬ 
nentially  stable  equilibrium  point  [9].  We  will  relax  this  condition  to  allow 
for  systems  with  multiple  invariant  sets.  Consider  a  switched  system  with 
C 1  functions  Vp  :  Mn  — >■  M  bounded  on  every  bounded  subset  of  Mn  such  that 


y(x)  <  0,  Vx  G  Rn, 
and  dp(x)  — >■  oo  as  x||  — >■  oo.  Denote 

(17) 

Ep  =  jx  G  Mn|l/(x)  =  o}  . 

(18) 

Our  problem  is  as  follows.  Is  there  a  dwell  time  condition  that  suffices  for 
ultimate  boundedness? 

Two  challenges  are  immediately  apparent.  First,  if  V  vanishes  outside 
of  M’p(k),  there  is  not  a  strictly  positive  Lyapunov  decay  rate  outside  of 
J\fp  (/%),  so  it  is  unclear  how  long  to  wait  between  switches.  Second,  even  in 
the  absence  of  switching,  the  active  system  may  never  enter  A fp  (k)  (Ep  may 
not  be  contained  in  Np(n)]  see  Figure  4).  The  following  three  subsections 
describe  two  methods  for  overcoming  these  challenges. 
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3.3.  Intermediate  Solution:  Expand  Entry  Neighborhood 

The  simplest  idea  is  to  realize  that  J\fp  (k)  grows  in  size  as  we  increase 
k.  Assuming  that  Ep  is  bounded,  we  can  pick  n  large  enough  so  that 
Ep  C  Mp  (k).  Even  so,  not  all  Lyapunov  functions  satisfying  (17)  will  decay 
exponentially  on  Mn  \  Nfp  (k).  For  example,  consider  a  single  one-dimensional 
subsystem  x  =  —  arctan(.r)  with  V  =  x2.  Nevertheless,  the  following  lemma 
is  useful: 


Lemma  2.  Consider  a  switched  system  with  Lyapunov  functions  Vp  satisfy¬ 
ing  Vp  <  0  with  Vp(x)  — >•  oo  as  ||x||  — >  oo  and  Ep  bounded.  Then,  for  sublevel 
Lyapunov  sets  Nfp(n)  such  that  Ep  C  Int{Nfp{n))  (where  Int  denotes  the  in¬ 
terior)  and  any  e  E  IZ,  there  exists  a  different  set  of  continuous,  radially 
unbounded  Lyapunov  functions  Vp  such  that, 

1.  Vp(x)  =  0  for  x  G  Nfp,  and  Vp(x)  >  0  for  xeR"  \  Nfp. 

2.  Vp(x(t  +  to))  <  e~£tVp(x(to))  for  all  x  G  Mn  \  Nfp  and  t ,  to  E  M. 

Proof.  Bhatia  [11]  constructed  a  unique  continuous  function  s(x)  on  Mn  \  Nfp 
such  that  s(x(t  +  to))  =  s(x(t0))  —  t  and  s(x)  — »  0  as  x  — >  Nfp.  Then,  we  can 
select  any  constant  £  and  set 


0  for  x  G  Nfp 
e£s(x)  for  x  G  Mn  \  Nfp. 


Hence,  on  Mn  \  Nfp, 

Vp(x{t  +  t0))  =  e£s(x(t+to))  =  eeWx(*o)M)  =  e-£tVp(x(to)). 


(19) 


(20) 


□ 


This  means  that  for  suitably  large  n,  there  exist  Vp  that  decay  exponen¬ 
tially  outside  Nfp  (n).  A  simple  constant  shift  can  patch  the  original  Lyapunov 
function  on  Nfp  (k)  with  the  construction  of  Lemma  2  outside  Nfp  (n)  at  dNfp 
(the  boundary  of  Nfp) ,  and  the  mere  continuity  of  Vp  outside  of  Nfp  does  not 
harm  any  essential  parts  of  the  proof.  (Note  that  performing  a  constant  shift 
outside  Nfp  (k)  will  scale  e  by  ap,  but  we  could  just  perform  the  construction 
again  with  a  larger  e  to  correct  for  this).  However,  relying  on  the  construc¬ 
tion  of  Lemma  2  may  not  allow  for  (10)  to  hold,  and  we  must  assume  that 
we  can  find  a  set  of  Lyapunov  functions  which  satisfy  both  the  exponential 
decay  property  and  the  //  property,  (10).  It  is  left  as  an  open  problem  to 
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determine  if  it  is  generally  possible  to  construct  a  set  of  Lyapunov  functions 
that  satisfy  both  conditions  for  any  switched  system. 

Now,  the  conditions  required  for  Theorem  1  hold  (noting  that  is  is  not 
essential  to  have  exponential  Lyapunov  decay  within  Mp). 

Remark  2.  It  is  useful  to  notice  that  our  set  definitions  do  not  require  a 
single  value  of  k  for  all  subsystems.  If  our  subsystems  have  qualitatively  dif¬ 
ferent  invariant  sets,  it  may  be  very  damaging  to  require  a  single  k.  Instead, 
we  can  choose  a  set,  n  =  {/%,},  of  values  and  restate  all  our  set  definitions  as 
functions  of  K. 

Making  np  very  large  is  problematic  in  two  ways.  First,  the  purpose  of 
defining  a  large  neighborhood  is  to  neglect  the  troubling  areas  (by  this,  we 
mean  that  we  do  not  know  much  about  the  trajectory  inside  A/”(k.)) .  However, 
with  large  kp,  we  may  be  cutting  out  substantial  portions  of  our  state  space. 
Secondly,  it  may  lead  to  larger  dwell  time  or  larger  F4(k).  In  the  next  two 
subsections,  we  describe  a  method  to  more  tightly  tailor  our  strategy. 

3-4-  First  Main  Result:  Tightly  Tailored  Entry  Set 

Since  the  troubling  regions  are  just  those  where  Vp  <  0  is  closest  to 
zero,  we  will  now  define  a  smaller  set  containing  these  regions.  Choose  a  set 
S  =  {rip! ,  ()p,2 .  ...,hPmax}  with  dPi  >  0,  and  define  the  set 

Up  {8P)  =  {xeK":yp  (x)  >  - 8P }  ,  (21) 

so  that  Ep  C  Ftp.  As  usual,  we  can  also  define  FL(S)  =  [J  FLP.  Figure  4 

pev 

provides  a  one-dimensional  example  to  help  visualize  how  these  sets  are  con¬ 
structed.  Similar  to  before,  not  all  Vp  satisfying  (17)  decay  exponentially 
outside  Fip  ( 5P ).  Define 


7p  (<$p)  =  maxDp(x).  (22) 

Furthermore,  set 


Cp  (Sp)  =  {xeRn:Vp  (x)  <  lp  (5P)}  .  (23) 

Since  Cp  is  a  sublevel  Lyapunov  set,  it  is  compact.  From  compactness,  Vp 
attains  a  minimum  on  Cp,  while  Vp  attains  a  finite  maximum.  Thus,  on  Cp  \ 
Ftp,  an  exponential  decay  rate  can  be  computed,  while  we  can  use  Lemma  2 
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Figure  4:  Qualitative  example  of  tighter  tailoring.  The  subsystem  has  two  stable  equilibria 
and  a  single  unstable  equilibrium.  Notice  that  J\fp  C  %p. 


outside  Cp.  Again,  if  necessary,  a  simple  constant  shift  can  patch  the  two 
functions  together  at  dCp. 

While  the  construction  in  Lemma  2  only  gives  s(x)  continuous  in  the 
multi-dimensional  sense,  it  is  clearly  directionally  differentiable  along  the 
subsystem  trajectories.  Thus,  writing  V  is  sound  notation.  Putting  it  all 
together,  we  can  compute  e  such  that 

Up  (x)  <  -eVp  (x)  ,Vxer\Hp  (5P)  ,  Vp  e  V.  (24) 

Unfortunately,  FLP  may  be  disconnected,  and  it  is  not  necessarily  invariant 
even  in  the  absence  of  switching.  We  will  engage  these  problems  directly  in 
Sec.  3.5,  but  for  now,  we  can  proceed  directly  to  a  simple  theorem  demon¬ 
strating  the  usefulness  of  embedding  A f  inside  FL. 

Theorem  3.  Consider  a  family  of  systems  defined  by  (2),  each  with  a  ra¬ 
dially  unbounded  Lyapunov-like  function  that  satisfies  (17).  Assume  FLP,  Ep 
bounded  and  Qp  C  Ftp(S).  Compute  s(S)  so  that  (24)  is  satisfied  and  assume 

Y^<riS),\/q,reV,  Vx  <E  Mn  \  "H  (<5)  (25) 

vq  lxl 

holds  for  finite  p(S).  Furthermore,  compute  np>  0  such  that  Qp  C  J\fp(n)  C 
HP(S).  Then,  for  every  switching  signal  with  dwell  time  r  >  los^^,  there 
exists  a  time  T  such  that  x(T“)  G  FL^t-)  (d). 
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Proof.  Notice  that  all  of  the  necessary  assumptions  are  valid  on  Mn\'Ha(t)(d) ■ 
By  taking  suitably  large  T,  we  can  make  14(t-)(x(T-))  arbitrarily  small. 
Therefore,  either  x(T_)  e  A fa(r~  )(«)  £  ^a{T~)  (S)  or  the  trajectory  enters 
PLa(t)  (3)  somewhere  else  before  that  time.  □ 

Remark  3.  In  [9],  k  was  a  single  tuning  parameter.  In  Section  3.3,  k  =  {kp} 
was  introduced  as  a  possible  set  of  tuning  parameters.  Now,  per  Theorem  3, 
the  set  of  tuning  parameters  is  5  =  {dp|,  and  n  is  computed  as  a  consequence 
of  our  selection  of  5.  Since  nearly  every  parameter /set  which  follows  is  de¬ 
pendent  on  8,  we  will  often  omit  explicit  dependence  in  favor  of  readability. 

3.5.  Second  Main  Result:  No-Escape  Set 

This  section  assumes  the  trajectory  has  entered  PLo-(t){S)  at  some  time 
and  proceeds  to  build  the  relevant  no-escape  set,  which  will  be  denoted  Aip. 

The  primary  problem  is  that  PLp  is  not  necessarily  an  invariant  set  even 
for  periods  of  time  when  a(t)  constant  is  constant  (i.e.,  no  switching).  For 
example,  a  subsystem  may  contain  a  locally  unstable  equilibrium.  With  5P 
small,  PL  is  certainly  not  an  invariant  set.  One  way  this  can  be  problematic 
is  that  we  can  get  a  free  extra  switch. 

Example  Consider  two  one-dimensional  subsystems,  iq  =  —x\  and  x2  = 
x2  —  x\,  with  V2  =  (x  —  l)2(x  +  l)2  and  small  <52-  Assume  a  nonzero  initial 
condition  with  the  first  subsystem  being  active.  The  trajectory  can  become 
arbitrarily  close  to  zero  before  switching  to  the  second  subsystem.  While  the 
second  subsystem  is  active  (. aft )  is  constant),  the  trajectory  can  clearly  leave 
PL>-  Furthermore,  since  the  second  subsystem  started  arbitrarily  close  to  the 
origin,  it  can  take  an  arbitrarily  long  time  to  exit  PL2-  Thus,  no  finite  dwell 
time  can  prevent  at  least  one  switch  from  being  possible  outside  of  PL. 

There  are  two  ways  to  compute  a  dwell  time  and  an  associated  spa¬ 
tial  bound,  but  we  need  a  few  more  definitions  first.  Set  the  usual  7(d)  = 

max  7 p(5p)  and  C  (S)  =  |J  Cp.  Compute 
pev  peP 

4  (S)  =  max  Fp  (x) ,  (26) 

xe£(<5) 

and  set  £(<5)  =  max  £P(S). 

pev 
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Theorem  4.  Consider  a  family  of  systems  defined  by  (2),  each  with  a  ra¬ 
dially  unbounded  Lyapunov-like  function  that  satisfies  (17).  Assume  Tip,  Ep 
bounded  and  Qp  C  TLP(5).  Compute  s(S)  so  that  (24)  holds  and  assume 


YL^<p(S),Wq,reV  Vx  £  E"  \  H  (S)  (27) 

Vq  fxJ 

holds  for  finite  11(d).  Furthermore,  compute  kp  >  0  such  that  Qp  C  Np(n)  C 
HP(S).  Set  k  =  min  np.  Then,  for  every  switching  signal  with  dwell  time 


T  > 


for  every  t  such  that  x  e  TLa(t),  x(f)  G  A4  (S) 


(28) 

U  Aip  for  all  t  >t,  where 

pev 


Mp(6)  =  {xeRn:Vp(x)<t;p}.  (29) 


Proof.  Consider  the  following  sequence  of  times.  Assume  a  switch  occurs  or 
the  system  is  started  at  to-  Assume  further  that  there  is  a  time  to  <  tenter 
at  which  the  trajectory  enters  'Hay+^(8).  There  may  or  may  not  be  a  time 
at  which  the  trajectory  exits  Tia^t+^(5),  which  we  will  denote  tenter  <  texit -  If 
there  is,  fexit  -  Center  could  possibly  be  arbitrarily  large.  Therefore,  the  second 
switching  time  t\  may  be  shortly  after  texit  while  the  trajectory  is  outside  of 
Ft^+AS)  regardless  of  the  dwell  time. 

The  proof  will  proceed  in  two  parts.  First,  we  will  show  that  the  trajectory 
does  not  leave  M.  before  a  third  switching  time,  called  t2.  Then,  we  will 
compute  a  dwell  time  so  that  the  trajectory  must  re-enter  Tia{t)  before  t2. 
By  induction,  the  trajectory  will  never  leave  M.. 


(i)  We  will  show  that  the  trajectory  never  leaves  M.  on  [tenter,  ti],  then 
on  [ti,t2].  Since  x(fenter)  G  FLay0)  C  Ma(to)  and  Ma(to)  is  a  sublevel 
Lyapunov  set,  the  trajectory  cannot  leave  M.a(t0)  on  [tenter  -  T]  since  no 
switching  occurs  on  this  interval.  Next,  notice 

xfa)  g  Ca{tT)  cCc  Ma(t+y  (30) 

Thus,  the  trajectory  cannot  leave  A4a(t+j  on  [ti,t2]. 

(ii)  We  will  show  that  the  trajectory  must  re-enter  FL^t)  before  t2.  Define 
A  =  t2-t1>T,  giving 
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The  dwell  time  given  by  (28)  gives  that  V^f-)(x)(i2  )  <  /%.  Thus,  the 
trajectory  enters  A /"CT(t+)  or  reenters  'Hcr(t)  elsewhere  before  time  £2 • 

□ 

Corollary  5.  Taken  together,  Theorems  3  and  4  provide  a  dwell  time  which 
guarantees  finite  time  entry  into  PLa(t)  and  the  invariance  of  A4  thereafter. 
Given  a  particular  problem,  if  we  want  to  apply  both  theorems,  we  should  set 

logmax(/i,|) 

T  >  - . 

£ 

Section  4  will  show  examples  that  have  ^  >  /i,  but  it  is  trivial  to  construct 
cases  with  very  large  //.  Next,  we  construct  some  lower  bounds  for  //  which 
are  helpful  for  determining  that  the  bounds  in  these  theorems  are  tighter 
than  some  other  possible  bounds.  Define 


VP(S)  =  maxVp(x), 

and  the  following  proposition  will  be  useful. 


(32) 


Proposition  6.  Given  the  preceding  set  definitions,  py  >  >  rjp  Vp. 

Proof.  Since  PL,,  C  Cp,  PL  C  £,  and  (p  >  rjp. 

Without  loss  of  generality,  consider  x  6  C  such  that  Vj  (x)  =  There 
are  three  cases: 

Case  1.  Suppose  x  G  Hi.  Then,  £i  =71.  Since  7  >  71  and  /1  >  1,  /17  >  £1. 

Case  2.  Suppose  x  G  TL  \  PL±.  Then,  £1  =  rp.  Since  £1  is  a  maximum  of  V j 
over  £  D  PL  and  H  is  open,  Vj  (x)  must  be  a  local  maximum.  But,  if  it  is  a 
local  maximum,  VVj(x)  =  0  and  Vj(x)  =  0.  Thus,  x  G  PL\,  and  the  problem 
is  reduced  to  Case  1. 

Case  3.  Suppose  x  G  £\PL.  Then,  (27)  holds,  and 


P  > 


U(x) 

U(x) 


Vp. 


(33) 


Then,  for  some  p,  py  >  pyp  >  pVj,(x)  >  Vj(x)  =  r/i  =  £1. 


□ 
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Notice  that  in  the  ideal  case,  with  each  subsystem  consisting  only  of  a 
single  exponentially  stable  equilibrium,  7  =  k,  Cp  =  Hf/  =  Afp,  and  £  =  a  can 

be  attained.  In  the  case  they  are  attained,  r  =  logma*(^’K)  _  logmax(y^)  _ 
,  (29)  collapses  to  (7),  and  we  can  consider  this  a  true  generalization  of 
the  result  in  [9].  In  fact,  one  can  notice  that  both  examples  in  [9]  satisfy  the 
equality  jin  =  £. 


Remark  4.  Using  similar  techniques  as  those  in  Theorem  4,  one  can  find 
a  bound  for  r  and  A4  that  looks  like  min(max(//7.  rjp),  £p).  This  expression 
reduces  to  given  Proposition  6. 

Remark  5.  We  chose  to  define  the  dwell  time  as  a  constant  across  all  subsys¬ 
tems,  but  allowed  the  spatial  calculation  to  proceed  with  regard  to  individual 
subsystems.  This  is  consistent  with  the  idea  that  applications  will  utilize  a 
single,  subsystem  invariant,  dwell  time.  Since  only  the  final  union  matters 
for  the  no-escape  set,  it  is  of  no  importance  that  the  individual  pieces  were 
built  with  different  constants  determining  sublevel  sets. 


4.  Examples 

We  provide  two  numerical  examples  to  illustrate  the  main  results  of  the 
paper.  In  particular,  we  want  to  demonstrate  that  tight  tailoring  can  provide 
better  results  than  expanding  the  entry  neighborhood. 

4-1-  One- Dimensional  Example 

Consider  a  family  of  systems  defined  by 

fi(x)  =  -x{x  -  l)2,  f2(x)  =  -(x  -  1).  (34) 

Select  Vi(x)  =  x 2  and  V2(x)  —  (x  —  l)2  and  5  =  [0.05;  0.1].  Then, 
Hi  =  [-0.139,0.197]  U  [0.803,1.139],  U2  =  [0.776,1.224],  //  =  67.28,  and 
£  =  1.29.  We  can  set  n  =  [0.0193;  0.050],  giving  Af\  =  [—0.139,0.139]  C  H± 
and  J\fr2  =  [0.776,1.224]  C  7^2-  Then,  we  can  compute  7  =  [1.297;  0.050], 
£1  =  [-1.139, 1.139],  £2  =  [0.776, 1.224],  and  £  =  [1.497; 4.575],  Finally,  this 
results  in  r  =  4.240  and  M.  =  [—1.224,  3.139].  The  size  of  M.  is  4.362. 

We  compare  the  large- k  method  with  two  metrics.  First,  we  select  n 
so  that  r  is  the  same  as  the  tight-tailoring  method  (we  make  the  methods 
temporally-equivalent)  and  compare  the  size  of  Ad.  Second,  we  select  n  so 
that  the  size  of  A4  is  the  same  as  the  tight-tailoring  method  (we  make  the 
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Figure  5:  The  sets  E  (black)  and  T-L  (red)  for  each  subsystem. 


methods  spatially-equivalent)  and  compare  the  required  dwell  time.  In  order 
to  make  r  =  4.240,  we  must  set  k  =  [2.271;  0.257].  This  results  in  n  =  8.837, 
£  =  0.514,  and  M.  =  [—1.507,3.507],  which  has  size  5.014.  For  the  second 
case,  we  can  make  A4  =  [—1.181,3.181]  by  setting  k  =  [1.395;  0.0328].  The 
result  is  that  /i  =  42.48,  £  =  0.0657,  and  r  =  57.08.  In  each  case,  tight 
tailoring  shows  better  performance. 


4-2.  Two-Dimensional  Limit  Cycle  Example 
Consider  a  family  of  systems  defined  by 


fi 


(x2  +  y2  —  1) 

1 


-1 

~(x2  +  y2 


f2 


(x2  +  y2  + 1) 

1 


-1 

-  (x2  +  y2 


f3 


(35) 


This  is  a  simplified  version  of  (1),  a  central  pattern  generator  used  for  control 
of  swimming  [12]  or  flying  robots  [2].  Subsystem  1  has  an  unstable  equilib¬ 
rium  at  the  origin  and  a  stable  limit  cycle  of  unit  radius  centered  around  the 
origin.  Subsystem  2  has  a  stable  equilibrium  at  the  origin.  Subsystem  3  has  a 
stable  equilibrium  at  [  2  0  ]  .  Take  the  corresponding  Lyapunov  functions 
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Figure  6:  The  sets  Q  (black)  and  J\f  (blue)  for  each  subsystem.  The  example  is  suitably 
well-behaved  so  that  E  —  Q  and  7~L  =  J\T  except  for  the  unstable  origin  and  its  surrounding 
neighborhood  when  p  =  1. 

to  be  Vi(x)  =  (x2  +y2  —  l)2,  V^(x)  =  (x2  +  y2)2,  and  V^x)  =  (P  —  2)2  +  y2)2, 
respectively.  If  we  choose  6  =  [0.2, 0.1, 0.1],  then  e  =  0.22,  /i  =  2215, 
k  =  [0.042,  0.022,  0.025],  7  =  [1.0,  0.22,  0.25],  £  =  [23,  33, 136],  and  r  =  38.9. 
Figures  5,  6,  and  7  show  the  resulting  sets.  In  fact,  many  of  the  result¬ 
ing  parameters  (chief  among  them  dwell  time  and  size  of  A4)  are  determined 
through  nontrivial  relations.  Optimization  and  tradeoffs  between  spatial  per¬ 
formance  and  temporal  performance  is  problem-specific. 

5.  Conclusion 

We  derived  a  stability  result  for  switched  systems  which  are  constructed 
from  subsystems  which  possibly  contain  multiple  invariant  sets.  This  amounted 
to  a  generalization  and  refinement  of  the  argument  presented  in  [9]  and  is 
in  the  spirit  of  dwell  time  methods  for  switched  systems.  This  result  can  be 
applied  to  a  larger  class  of  dynamical  systems  than  those  in  [9],  including 
those  which  contain  bifurcations.  We  provide  a  set  of  tuning  parameters,  S, 
which  cause  interrelated  and  problem-specific  effects  on  performance. 

Future  studies  could  be  focused  on  determining  an  explicit  relation  be¬ 
tween  tuning  parameters  and  performance  for  a  subclass  of  systems  in  or¬ 
der  to  better  facilitate  optimization  or  investigating  a  continuous  family  of 
switched  systems  near  a  bifurcation  point. 
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Figure  7:  The  sets  Ep  (black),  £  (red),  and  A4  (blue).  The  skinny  black  line  is  an  example 
trajectory. 
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